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Introduction 


This is a modelling unit. Its subject is one which is not only of interest to many 
people, but also lends itself to modelling particularly well, in that even the simpler 
models can be made to give useful answers to practical questions. It will also give 
you the opportunity of using, in a practical context, some of the mathematical 
techniques discussed earlier in the course. If you have studied TM281 you will find 
some of the topics familiar. 


Study guide 


This unit is intended to be read as it is written—the sections following each other 
in numerical order. The first four sections are normal teaching sections; Section 5 
contains no new material, but consists of numerical problems somewhat more 
difficult and comprehensive than the exercises in the first four sections. You should 
attempt as many of these problems as you can without looking at the solutions. 
When you have done all you can (or all you have time for), read the solutions 
carefully—you should find them instructive. 


There is no tape associated with this unit but there is a television programme. The 
exact stage of your study of this unit at which you watch it is not too important, 
except that you should have read at least Section 1 first. The ideal point to watch 
it is when you reach Section 3. 
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1 The nature of heat transfer 


1.1 What is heat transfer? 


This unit is about the transmission of energy from one place to another. For an 
example you might look at the way in which the energy produced by the sun is 
available to plants on earth, or the way in which some of the energy from a gas 
fire, intended to heat a room, escapes through the walls to heat the air outside. 
This kind of redistribution of energy is conventionally referred to as heat transfer. 


Actually, this term is something of a misnomer. It leads to such misleading 
statements as ‘heat flows from one body to another’—as though we were dealing 
with a material substance; a kind of fluid perhaps. A more accurate, though 
lengthier, title for this unit would be Energy transfer by heating, where ‘heating’ is a 
particular mechanism for acquiring or disposing of energy. 


The essential feature of this mechanism is the existence of a temperature difference, 
combined with the rule that the net transfer of energy takes place from regions of 
higher temperature to regions of lower temperature. For example, the fact that 
energy is transferred from a heated room to the cooler atmosphere outside the 
house, through the walls, is only too clear to the householder who has to pay for 
the fuel which keeps the heater going. Again, the energy reaches the air in the 
room via, say, a central heating radiator which is itself at a temperature above that 
of the room. 


These are just two examples of the kind of process we shall try to model in this 
unit, and if now and again I use the term heat transfer, it is because it has become 
generally accepted. It will not be allowed to introduce any unnecessary difficulties. 


Heating is not the only way in which energy can be transferred. If you wind up a 
watch you are transferring energy from your muscles to the watch spring. A 
temperature difference is irrelevant to this operation (known as transferring energy 
by working), and it will not therefore be dealt with in this unit. However, the fact 
that the discussion will be restricted to energy transfers which rely on temperature 
differences for their driving force does not mean that our examples will lack 
variety. There are in fact three distinct modes of heat transfer, and I shall have 
something to say about each of them. 


I shall develop some specific models which are widely used, sometimes to predict 
the rate at which energy is being transferred (or the amount of energy transferred 
in a given time), and sometimes to predict the way the temperature will vary inside 
a piece of material. First of all, however, I must say something about the three 
modes of transferring energy by heating. 


1.2. The three modes of heat transfer 


The three modes of heat transfer are called conduction, radiation and convection. In 
practice they often occur in combination, but for our present purpose it is 
convenient to consider them one at a time. 


In conduction, energy is transferred from one part of a piece of material to another 
part because the faster movement of atoms or molecules or electrons in the hotter 
part (i.e. the part at higher temperature) excites the atoms or molecules or 
electrons in the cooler part so that they move faster. In this way energy is passed 
on along the body from the parts at higher temperature to those at lower 
temperature. There is no macroscopic (i.e. large-scale) bodily movement of any 
part of the material. I have already mentioned one example of conduction: the loss 
of energy through the walls of a heated room. The inside face of the wall is at a 
higher temperature than the outside face and, as a result, energy passes through 
the wall by conduction. The same process can take place between two different 
pieces of material in contact; energy will pass from the hotter to the cooler. 


The second mode of heat transfer is radiation. A familiar example of radiation is 
the way that energy reaches us from the sun. Unlike conduction, radiation does 
not require the presence of any material object. Radiation travels through empty 
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space. The mechanism by which the energy travels is called electromagnetic 
radiation. This course is not the place for a detailed discussion of electromagnetic 
radiation, but it is worth mentioning that radio and television broadcasts, light 
and X-rays are all examples of it. But radiation does not take place only in empty 
space. Every piece of material continuously dissipates part of its energy by 
radiation and gains energy by radiation from the other objects around it. For a 
given body the net rate at which it gains or loses energy by radiation depends on 
its temperature and that of its surroundings, as well as the nature of its surface. 
The higher the temperature of the body the greater the amount of energy it 
radiates in a given time. A highly polished surface is a poor radiator and also a 
poor absorber of radiation. This is why in a Thermos flask the surfaces which face 
each other across the vacuum space are highly polished reflectors. 


The third method of energy transfer, convection, is associated with fluids (i.e. 
liquids and gases). Energy is transferred by virtue of the movement of the fluid. A 
saucepan of water on a gas ring is a suitable example. Energy from the burning 
gas is conducted through the bottom of the saucepan and transmitted to the parts 
of the water nearest the bottom of the pan, mainly by conduction. This raises the 
temperature of that part of the water and consequently makes it less dense. The 
hotter water then rises and mixes with the main body of the water, transmitting 
energy to it, again mainly by conduction. It is replaced near the bottom of the pan 
by cooler water from above and the process continues. As a result, there are 
convection currents in the water. If they are caused purely by the kind of variation 
in density with temperature that I described in the case of the water in the 
saucepan, the process is known as free convection. If the fluid movement, and 
therefore the rate of heat transfer, were to be increased by the use of a stirrer or a 
pump, then it would be forced convection. 


1.3. SI units concerned with heat transfer 


Before we begin to discuss the modelling of the heat transfer processes, | want to 
introduce an important relationship between energy and temperature. It tells us 
the amount of energy that is required to produce a given change in the 
temperature of a fixed amount of a substance, and we shall need it later on in this 
unit. 


Think of a quantity of water being heated in an electric kettle. Electrical energy is 
being used up by the heating element (as the meter will testify). What happens to 
this energy? It goes to raise the temperature of the water and is therefore stored in 
the water as long as the water stays hot. While the water stays hot the energy 
originally supplied to it by the heating element is available for transfer elsewhere; 
to heat a cold bed by way of a hot water bottle, for example. The energy stored in 
a substance by virtue of its temperature is known as the internal energy of the 
substance; I shall denote it by the symbol E. This is another form of energy to add 
to those which you have already met: for example, the kinetic energy of a body (by Unit 4 
virtue of its velocity) or its potential energy (by virtue of its position in, say, a 
gravitational field). The change in the internal energy of a mass m of water in 
raising its temperature from 0, to @, 1s found to be 


E, —E, =mc(@, — 6;) (1) 


where E, is the internal energy when the temperature is 0, and E, that when it is 
0,. The quantity c (which may, for the purposes of this unit, be considered 
constant) depends on the nature of the substance being heated, and is called the 
specific energy capacity or sometimes the specific heat of the substance. 


Two points about Equation (1) are worth noting. The first is that this equation 
does not describe the process of heat transfer to the substance, but only the result 
of such a process. The second point is that this equation applies to any kind of 
substance—it does not have to be a liquid. It breaks down in complicated cases 
like mixtures of gases and liquids (e.g. when the water in the kettle begins to boil), 
but these cases are beyond the scope of this unit. For gases the relevant value of c 
is that which corresponds to heating at constant volume. 


-— Silvered surfaces 
(facing vacuum) 


Asbestos 
anti-vibration pad 


Vacuum seal 


This brings me to the units we are going to use. As you know, we use the SI See Unit 4 and the Handbook. 


system of units in which energy is measured in joules (J) and mass in kilogrammes 
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(kg). We also need temperature, which is measured in degrees Celsius (°C) 
(commonly called centigrade). What then are the units of specific heat? 
From Equation (1), we can put 
E, == 1 
eo m(8, — 8,) 
Substituting the units of the terms (E, — E,), m and (6, — @,), we find that the 


J ; 
units of c are ——~ or Jkg™ ‘°C * (in words, joules per kilogramme per degree 


ke °C 
Celsius). 
Approximate values of c in these units for a number of substances are shown in 
Table 1. 
Table 1 
Specific heat 
Substance Ate Ce The oblique stroke means 
‘expressed in units of”. 
water 4200 
alcohol 2 500 
copper 390 
glass 630 
air heated at 715 
hydrogen constant volume 10000 


There is one further quantity which we shall come across in the course of our 
work, and that is the rate at which energy is produced (e.g. by the heating element 
in the kettle). The unit of time in the SI system is the second, and so from what I 
have already said it is clear that the rate of energy production (or transmission or 
dissipation) is measured in joules per second (Js~*). This unit has a special name: 
it is the watt (W). Hence 1W = 1Js“'*. 


If, for example, the element in the kettle is rated at 3 kW, it will use electrical 
energy at the rate of 3000 joules per second, as long as it is switched on. (This rate 
of energy use is called the power of the heating element.) 


Here are a few exercises to enable you to practise using the quantities and units I 
have mentioned in this section. 


Exercise 1 


Work out the change in the internal energy of 2 kg of water when its temperature is raised 
from 20°C to 80°C. Use Table 1. 


Exercise 2 


The water in Exercise 1 was heated in an electric kettle with an element rated at 2 kW. 
Assuming that all the electrical energy is used only in heating the water, estimate how long 
it will take to raise the water temperature from 20°C to 80°C. 


Exercise 3 


In Exercise 2 it is assumed that the electrical energy is used only to heat the water. What 
else would it be used for, in practice? Is the answer to Exercise 2 an underestimate or an 
overestimate? 


[Solutions to Exercises 1-3 on p. 27] 
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Summary of Section 1 


Heat transfer means energy transfer by heating. This kind of energy transfer 
requires the existence of a temperature difference, and the energy is transferred 
from regions or bodies of higher temperature to those of lower temperature. There 
are three modes of heat transfer: conduction, radiation and convection. 


Part of the energy of a body depends on its temperature. This part is known as 
internal energy, and the change of internal energy due to a change in temperature 
(0, — 0,) is mc (0, — 0,), where: 


m = mass of body, 
c = specific energy capacity or specific heat. 


2 Steady-state conduction in one dimension 


2.0 Introduction 


It has long been known that when different parts of a piece of material are at 
different temperatures, energy is transferred from the regions of higher temperature 
to those of lower temperature. The external walls of a house are a case in point. If 
the inside of the house is warmer than the air outside, then the inside surface of 
the wall will be warmer than the outside surface and energy will be transferred out 
of the house through the wall by conduction. It will then be transferred from the 
outside of the wall to the air, mainly by convection. In this section I want to 
concentrate on thermal conduction through solid walls. 


The model I shall use depends on several simplifying assumptions: 


(1) The wall can be represented as a uniform flat slab. In consequence we can 
choose coordinates so that one face of the slab is in the plane x = 0 (see Figure 1) 
and the other is in the plane x = b, where b is the thickness of the slab. 


(2) The slab is in a steady state, which means that the temperature at any given 
point is independent of the time. 


(3) The temperature in the slab is uniform in each plane x = constant. This 
implies that the temperature in the slab at any given time depends only on x and 
not on y or z. When the temperature depends on only one space coordinate in this 
way, we say that the temperature distribution is one-dimensional. 


Assumptions (2) and (3) imply that the temperature in the slab, being 
independent of time t and the space coordinates y and z, depends only on the 
space coordinate x. Writing 0 for the temperature we can say that 6 is a function 
of x only, or in symbols 


6 = 0(x) (0<x <b). 


The function 6(x) gives a complete description of the temperature distribution in 
the slab. The question now is how to find this distribution and how to find the 
rate of energy transfer across the slab. 


2.1 Fourier’s law 


The model that we shall use for steady-state conduction in one dimension is based 
on theoretical work done by Joseph Fourier in the early nineteenth century which 
has since been amply confirmed by practical experience. Fourier’s law gives 
quantitative expression to the idea that heat conduction transfers energy from 
hotter to cooler places, by postulating that the rate of energy transfer is 
determined by how rapidly the temperature varies with position. 


The law gives the rate of energy transfer by conduction across any small area, 
called an area element, within any block of material. For the one-dimensional 
problems considered in this unit, we need only consider area elements which are 
parallel to one of the coordinate planes. If the small area is parallel to one of the 
coordinate planes y = 0 or z = 0 (see Figure 2), then since the temperature does 


Figure 1 


Plane x = 0 


Plane x = b 


Here were are using three- 
dimensional Cartesian 
coordinates (x, y,z). They are 
explained in more detail in 


Unit 14. 


Small area 
parallel to 


v=o a 


Figure 2 


Small area 
parallel to 
x=0 


Small area 
parallel to 
z=0 
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not vary in the y and z directions, we may expect that no energy is transferred 
across it. On the other hand, for an area element parallel to the plane x = 0 (ie. 
parallel to one of the faces of our slab), there will be transfer of energy. 


We may expect the rate of heat transfer (q) across this area element to be 
proportional to its area, which may be denoted by a. The rate of heat transfer will 
also depend on how rapidly the temperature changes with position, that is to say 


d ; 
= which is called the temperature gradient. Fourier made the natural assumption 


ek do 
that q is simply proportional to ae Combining these two proportionality 


relations, we are led to write 


oc a—. 
q ie 


Before writing this as an equation, we should consider the sign of the constant of 


proportionality. If ma is positive, the temperature :increases as x increases; and 
since heat flows from hot places to cold the heat transfer in this case will be from 
large values of x to small, i.e. in the negative x direction. The natural sign 


convention for g, however, is that a positive value of g corresponds to flow in the 


re, dé es 
positive x direction. So — and q have opposite signs. To allow for this, we write 


dx 
Fourier’s law as an equation in the following form: 
dé 
uae il 1 

a= ke _ (1) 
or 

ee ve 

a ds 


where k is a positive constant of proportionality. 


The constant of proportionality « in Equation (1) is called the thermal conductivity 
of the material through which conduction takes place, and different materials can 
have very different values of x. Although Equation (1) assumes x to be a constant, 
many materials are found in practice to have a value of x which varies with 
temperature. For the purposes of this unit we can neglect this effect. Before we 
look at representative values of x, let us just be clear about the SI units of thermal 
conductivity. 


Question What are the SI units of thermal conductivity? 
Answer From Equation (1) we have 


a 
ad6/dx 
Now the SI units of g are watts or W, the SI units of a are (metres)* or m7, and 


° ° 


or —. The negative sign is irrelevant to the question 
m 


dé 
the SI units of — are 
dx metre 


of units. 
ae W 
a Cin. m°C 


or Wan 3c -. 


So the units of x are 


Table 2 shows some typical values of « for various solid materials, all of them 
fairly common. Notice that the ratio of the largest value (for copper) to the 
smallest value (for cork board) is nearly 9000, and even for the two metals (copper 
and mild steel) the ratio is about 7. 


10 MST204 12.2 


Table 2 
Thermal conductivity 

Material wwa“Cc™ 
copper 380 
mild steel — 54 
concrete 1.4 
glass 1.0 
brick 0.7 
wood (teak) 0.17 

_ asbestos board 0.16 
cork board 0.043 


2.2 Finding the temperature distribution 


Before we can use Fourier’s law to determine the temperature distribution 6(x), we 
need to know how gq depends on x. To find this out, we now consider not an area 
element but a volume element, that is a small region inside the slab, bounded by 
faces parallel to the coordinate planes (see Figure 3). 


We can apply the input-output principle (Unit 3) 
accumulation = input — output 


to this volume element. First of all, since we are assuming a steady state, the 
temperature of the volume element does not change and so its energy does not 
change. That is, there is no accumulation, so input and output must be equal. Our 
volume element has six faces, but there is no heat transfer across the four faces 
that are parallel to the coordinate planes y = 0 or z = 0 because the temperature 
does not vary in the y or z directions. The only heat transfer is across the two 
faces that are parallel to the x = 0 plane, which are shaded in Figure 3, and since 
input and output are equal, the heat transfer rates across these two faces must be 
equal. That is, for a steady state the values of g for the two opposite faces of a 
volume element are equal. With the geometry we are considering here, the areas of 


those faces are also equal, so 4 is the same at the two ends of the volume element. 
a 


dé 
It follows, by Fourier’s law, that Ke is the same at the two ends of the volume 
x 


element. 


The above argument applies to any volume element in the slab; for example, we 
could take one end of the volume element in the face x = 0 and the other end in 


dé 
some arbitrary plane x = x,. Then the argument tells us that Ke has the same 
x 


value at x = x, as it does at x = 0. This is true for all x,, and so we have shown 


da 
that —(x) is a constant. 
dx 


Let us call this constant B; we then have 
dO 
a 
This differential equation can be solved by direct integration to give 


06=Bx+C 


B. 


where C is another constant. This means that the graph of 0 against x is a straight 
line. 
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To determine the constants B and C we need further information. For example, 
suppose we know the temperatures of the two faces of the slab to be 6, and 0,, 
say. Then we have 


0d, = 000) = C, 
0, = 0(b) = Bb+C 
where b is the thickness of the slab measured in the x direction (Figure 3). Solving 
for B and C gives 
te ae 
ae wee 
C= 0,. 


B 


Once we know the temperature distribution we can calculate the rate of heat 
transfer through the slab. By Fourier’s law the heat transfer per unit area is 


ed 
es ae 
= —KB 
0, —9, 
= —K . 


2.3 Conduction through a constant cross-sectional area 


We are now ready to model steady-state conduction through a wall. I shall assume 
that the inside of the house is warmer than the outside and that energy transfer 
takes place only from the inside face of the wall to the outside face; I shall neglect 
any energy that passes to the rest of the structure via the edges of the wall. I shall 
also assume that the temperatures inside and outside are uniform and that the 
wall is made of uniform material like the slab I mentioned earlier. What we want 
our model to tell us is the rate at which energy passes through the wall. 


In terms of the slab shown in Figure 3, the assumption that we can neglect any 
energy transfer through the edges of the wall implies that the temperature will not 
vary in the y or the z directions, and so we are justified in using the model 
described in Subsection 2.2. The total heat transfer rate will be the sum of the heat 
transfer rates of the elements, i.e. at any cross-section of the wall perpendicular to 
the x direction, 


d 
ee 3 [-a‘e] (where ) stands for ‘the sum of ...’) 


dé 
= —K7 da 


dé 
since kK and Fe are the same for all the elements. 


It follows that 64 Direction of 


heat transfer 
Ave eee ead 
Gotal — —~KAG_ 


dx 
where A is the area of the wall at right angles to the x direction, so that )' a= A. 


Figure 4 shows the temperature variation in the x direction. Here, 0, is the 
temperature of the inside face and 6, is the temperature of the outside face, so that 
G, >. 

0, ies 


As before, - = “a where b is the thickness of the wall, so that 


0, —90 
Qtotal = —KA lores | | 
b 


= KA a5. (2) Figure 4 
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Equation (2) is the most commonly used version of the steady-state conduction 
equation. The following exercises are based on it so that you can see that it 
supplies the information we set out to obtain at the beginning of this subsection. 


Exercise J 


The inside face of a brick wall 240mm thick has a temperature of 18°C, and its outside face 
has a temperature of 10°C. Estimate the rate of heat transfer per m? of wall surface. (You 
will need to use Table 2.) 


Exercise 2 


What is the temperature gradient for the wall in Exercise 1? What is the temperature at a 
section of the wall 100mm from the inner face? 


[Solutions to Exercises 1, 2 on p. 27 | 


2.4 Conduction through a changing cross-section 


I want to end this section by discussing the modelling of conduction through pipe 
walls—for example, central heating pipes which carry hot water to the radiators 
and are surrounded by cooler air. I shall again assume a steady state, as defined 
earlier, and also that heat transfer takes place radially (i.e. at each point in the « 

pipe wall, the direction of heat transfer points away from the centre of the pipe). 

The main difference between this case and those we met in earlier pages is that 

this time the cross-sectional area through which energy is transferred is not : 
constant. This becomes clear when we consider the cross-section of a pipe such as ere 
that shown in Figure 5. The pipe is of length |, and has inner radius r, and outer 
radius r, (so that the thickness of the pipe wall is r, — r,). 


Consider a cylindrical surface of radius x inside the pipe wall (so that r, < x <r,), 
as shown in Figures 6(i) and (ii). 


Its area is 2xxl. By symmetry, the temperature at all points of this surface will be 
the same. It follows that the steady-state temperature @ in the pipe wall will be a 

: dé ' 
function only of x and the temperature gradient will be ax If we consider a small 
part of the curved surface, it will behave just as if it were flat, and so the rate of 

: do 

heat transfer per unit area is, as in the case of the flat slab, equal to ae Hence 
if q is the steady-state heat transfer rate through the pipe wall and A is the area of 
the cylindrical surface, 


<2 
A dx 
or 
dé 
q = KA 


where A is not constant, but a function of x. In fact, as we have seen, A = 2zxl, 


and so 
a= (0) x Ors x (i) — 
dx 
Hence — l 
dé —q 


Or a 


dx 2nkxl' Figure 6 


In the steady-state, q is a constant. Thus the right-hand side of Equation (3) is a Unit 2, Subsection 3.1 
known function of x, and this equation can be solved by direct integration to give 
the general solution 


oe = 
ais let “ib 


q 
K 


— 57] OBe * +C (4) 


where C is a constant of integration. 
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Example 


Sketch the temperature variation along a radius of a pipe carrying hot liquid. Will 
the temperature halfway between the inside and outside surfaces be greater or 
smaller than it would be for a flat wall of the same material and thickness with the 
same inside and outside temperatures? 


Solution 


We use the notation of Figure 6(i). The inside temperature 6, will be the greatest. 
The gradient of the graph of 0 against x will be negative at all points but 
flattening out as x increases. This follows from Equation (3). The required graph is 
sketched in Figure 7 (full line). The graph for a flat wall is shown dotted. Clearly 
the half-way temperature is less for the pipe than for the flat wall. 


By substituting into Equation (4) the values appropriate to the inside and outside 
curved surfaces of the pipe, we get 


ee 
0, = Fad Oee™! +C 
and 
Ua ae =selog r+ C. 
: onal" * 


Subtracting the second from the first, we have 


q 
0, — 8, = 5 eee r, — log. r) 
=~" (log, r, — log, r) 
Ps eee . 


SO 
_ 2nkl(0, — 9.) (5) 
tog. (2/74) 


Exercise 3 


A glass pipe has an inside diameter of 50mm and an outside diameter of 80mm. The inside 
surface of the pipe is at a temperature of 90°C and the outside surface is at 50°C. Estimate 
the rate at which energy is conducted through the wall of the pipe per metre length of pipe, 
in the steady state. 


Exercise 4 


For the pipe specified in Exercise 3 what is the temperature in the pipe wall at a radius of 
32mm? 


[Solutions to Exercises 3 and 4 on p. 27 | 


Summary of Section 2 
Fourier’s law for steady-state conduction in one dimension is 


= —Ka— 
. dx 


where: g = heat transfer rate in the positive x direction, 
k = thermal conductivity, 
a = area at right angles to the direction of heat transfer, 
= temperature, 
x = distance measured in direction of temperature variation. 


dé 
ae is known as the temperature gradient. 
x 


Figure 7 


Remember that g is constant 
and has the same value at 
every radius. 
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For steady-state conduction through a uniform slab of constant cross-sectional 
area A, Fourier’s law reduces to 
(0, ae 94) 


KA 
ge b 


where (6, — 95) is the temperature drop across thickness b of the material. 


For steady-state conduction through a pipe of circular cross-section of internal 
radius r,, external radius r, and length |, Fourier’s law leads to 


_ 2nkl(0, — 8) 
loge (r/r) 
where (6, — 8,) is the temperature drop across the pipe wall. 


3 Television programme notes 


Now watch the television programme: ‘One-dimensional steady state heat transfer’. 


TV 12 


Read Section 3 after viewing the programme. (This section is intended to be read 
after watching the television programme, but you should read it and work through 
the numerical parts even if you missed the programme, because it contains new 
assessable material.) 


The main purpose of the television programme is to investigate the practical 
usefulness of the simple models of heat transfer which are derived in the unit, 
especially that for conduction. Dr Audrey Stuckes is studying the thermal 
insulation properties of building materials. She is looking at the performance of 
fibrous roof insulation materials for the Building Research Establishment and at 
the properties of polymeric materials in wall structures for the Polymer 
Engineering Directorate. The experiments comprise the careful measurement, 
under realistic conditions, of temperatures and heat transfer rates. Although 
conditions in practice are rather more complicated, it turns out that the heat 
transfer in buildings can be adequately modelled by one-dimensional steady-state 
models. Indeed, if the models are to be really useful to builders and architects, they 
must be simple. 


In order to remind you of the simple model for conduction and to introduce the 
measuring devices used in the work described so far, the programme includes a 
simulation of one-dimensional heat conduction along a metal bar which has a 
heater at one end and heat flow meter at the other. (We have used a simulation 
because we wanted to remind you of the fact that when the heater is first switched 
on the temperature at any point in the bar will change with time but will 
eventually settle down to a steady value, and the time required to do this in an 
actual experiment would have been excessive.) 


The main point, made with the help of animated diagrams, is to illustrate what is 
meant by ‘one-dimensional steady-state heat transfer by conduction’. 


Exercise 1 (revision) 


The metal bar mentioned above is 12mm in diameter and 200mm long. The steady-state 
temperature difference between its ends was shown as 10°C and the heat flow meter showed 
a reading of about 1.12 W. Estimate the thermal conductivity of the material of the bar. 


[Solution on p. 27] 


Although conduction is the most significant mode of heat transfer in building 
materials, convection and radiation also play their part. Yet the total heat transfer 
rate is, in practice, usually modelled by an even simpler expression than that for 
conduction. This is done by using a quantity called the U value, which is a kind of 
overall heat transfer coefficient. It can be, and is, used in cases where heat transfer 
occurs across more than one layer of material (as for example, in a cavity wall). 
The U value in such a case is worked out by dividing the heat transfer rate per 
unit area by the total temperature difference across the material in question, and it 
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can take account not only of conduction through solid materials but also of 
convection (as in the heat transfer between a wall and the air on either side of it) 
and radiation effects. Here is an example. 


Example 1 


The air temperatures on either side of a cavity wall are 21°C and 0°C respectively. 
It is found that the total heat transfer rate per square metre of wall is 30 W. What 
is the U value for the wall? 


Solution 
_ heat transfer rate per unit area 
total temperature difference 


= 3) = 145 Wm ?°C™'. 


U value 


I shall have more to say about U values in Section 4. For the moment it is enough 
to know that the U value is a measure of the ease with which heat transfer can 
take place and that at present, because of the need to conserve energy, there is a 
premium on methods of building construction which result in low U values. 


The programme also features demonstrations (not simulations) of convection and 
radiation. Convection is shown by the movement of particles in a liquid which is 
heated at one point. Radiation from an electric heater, via a reflector, is used to 
ignite a match. A very important fact about radiation is mentioned, namely that 
the rate of energy emission by radiation depends on the fourth power of the 
absolute temperature of the radiating body. This needs a little more explanation, as 
follows. 


Radiated energy is emitted and absorbed by any given body all the time. In this 
unit we deal only with cases where the absorbed radiation is very much less than 
the emitted radiation (for example, an electric lamp) and where we can say with 
little error that the energy emitted by radiation equals the net heat transfer rate, q. 
Now in such a case, 


q = y0s 
where yp is a constant and 6@, is the absolute temperature, defined as follows: 


where @ is temperature in degrees Celsius. This definition implies that the zero of 
the absolute temperature scale is at — 273°C. This point represents the lower limit 
of attainable temperature, and so the absolute scale of temperature is not based on 
the properties of a particular substance (in the way that 0°C refers to the freezing 
of water, for example). 


The SI unit for absolute temperature is the Kelvin, represented by the symbol K 
as in the statement ‘0° C is equivalent to +273 K’. Note: (i) that one writes ‘K’ 
and not ‘degrees K’; (ii) that the temperature intervals are the same on the Kelvin 
scale as on the °C scale. (For example, the temperature difference between 240 K 
and 230K is the same as that between 240°C and 230°C.) This is why the units 
of, for example, thermal conductivity are often quoted as Wm~! K~", which is 
exactly equivalent to Wm~‘4°C™~* since temperature differences are relevant here 
rather than particular values of temperature. 


It is also worth remembering that, as I said in Section 1, heat transfer by radiation 
can take place across empty space. 


Exercise 2 


Write down the absolute temperature, on the Kelvin scale, which corresponds to each of the 
following: 


(i) —250°C, Gi). 3°€. 
Exercise 3 


A vacuum lamp is rated at 40 W. Write down an expression for the heat transfer rate from 
the lamp as a function of the filament temperature 0°C and estimate the value of y if 
6 = 2000. 


[Solutions to Exercises 2 and 3 on p. 27] 
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Summary of Section 3 


The television programme includes demonstrations of conduction, convection and 
radiation, and of research work into the thermal insulation properties of building 
materials carried out for the Building Research Establishment and the Polymer 
Engineering Directorate. 


The heat transfer through walls and roofs is in practice often represented by 
simple, steady-state, one-dimensional models. In particular the U value is in 
common use as an index of how readily heat transfer can take place. The U value 
is defined as follows: 


heat transfer rate 


U value = ——____________., 
overall temperature difference 


New techniques and components are being developed to satisfy the demand for 
very low U values in buildings. 


The thermal energy emitted by radiation from a body depends on the fourth 
power of the absolute temperature of the body, i.e. 


thermal energy emitted by radiation = y6;, where: 


uu is a constant, 
0, = 273 + 6 = absolute temperature, 


0 = temperature in °C. 


4 Convection and insulation 


41 A simple model of convection 


So far, we have looked at conduction through solid walls without paying much 
attention to the way in which the energy is transmitted to and from the wall. This 
is what I want to consider now. In the case of the house wall, energy is transferred 
from the warm air in the room, through the wall, to the cool air outside. There are 
therefore two regions where there is heat transfer between a fluid (air) and a solid 
(the wall). In the case of the pipe, the energy is transferred from the hot water 
inside to the cool air outside, and again there are two fluid—solid boundaries 
(water—pipe and pipe—air) across which heat transfer takes place. In either case, 
since we are dealing with a fluid, convection will play a dominant part in the heat 
transfer process associated with the boundary. 


Considered in detail, convection turns out to be a rather complicated matter, not 
at all easy to analyse or to model. As you would expect from what I said in 
Section 1, fluid mechanics plays a very important part in any thorough 
consideration of convection. Fortunately, in spite of the fact that fluid mechanics is 
outside the scope of this course, it is possible to get some quite useful results based 
on a very simple description of the main features of convection, together with 
some of the practical experience that has been accumulated over the years. 


Figure 1 shows a solid—fluid boundary, and I shall assume that heat transfer is 
taking place from the solid (say a house wall) to the fluid (atmospheric air). As we 
know, this necessarily means that the temperature of the face of the wall, where it 
touches the air, must be higher than the air temperature at some distance from the 
house. The parts of the air near the wall are heated by conduction and radiation 
before they move away to join the main body of the air, being replaced by cooler 
air which in turn is heated. The result is that the air very near the wall has a 
temperature higher than that of the main body of the air, which is assumed to be 
at a uniform temperature. So the temperature drop between the wall and the main 
body of the air, which enables heat transfer to take place, is concentrated in a thin 
layer of the air quite near the solid boundary. This is sketched in Figure 1, where 
O, represents the temperature of the wall at its boundary and 0, represents the 
temperature of the main body of the air. 


Figure 1 


Fluid (air) 


0; 


Graph of temperature 
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The same general argument applies when the direction of the heat transfer is from 
the fluid to the solid; the temperature drop from fluid to solid is again 
concentrated in a thin layer of fluid immediately next to the solid boundary. In 
this layer the fluid is giving up some of its energy to the wall, mainly by 
conduction. Figure 2 shows part of the cross-section of a solid house wall, together 
with a sketch of the temperature variation from the inside air at temperature 0,, to 
the outside air at temperature 0,,,. We know how the temperature drop (0, — 0,) 
across the wall is related to the heat transfer rate—but what about (0, — 0,,,)? It 
is reasonable to assume that the nature of the solid surfaces and the nature of the 
fluid will have something to do with it. We also need to know whether the fluid 
near the wall is stationary or moving. (Inside the room the air is presumably 
stationary or at most moving very slowly; outside the house there may be strong 
winds.) 


The model for convection that is commonly used in practice by chemical 
engineers, architects and heating engineers is very simple. It is as follows: 


gq =hAAO@ (1) 
where A = area of the solid—fluid boundary, 
A@é = temperature difference between solid and fluid, 


q = rate of heat transfer, 
and h is a quantity called the convective heat transfer coefficient. 


Equation (1) looks similar to the equation we used to model conduction through a 
K é 
slab, with h substituted for b In fact, h is a different sort of coefficient because it 


takes into account not only the nature of the solid surface but also the nature of 
the adjacent fluid, as well as the relative velocity between fluid and solid surface. 
Question 

What are the SI units of the convective heat transfer coefficient h? 

Answer 


oo 
sae 


or Wm *°C:'. 


so that the units are 
rs glee @ 


Question 


In the case of the house wall in Figure 2, will the value of h be greater for the 
inside surface or the outside surface? 


Answer 


The greater movement of the air outside will speed up the process of convection 
and hence the heat transfer rate. Consequently h will be greater for the outside 
surface than the inside surface. 


The actual value of h to be used in Equation (1) is based on measurements of heat 
transfer rates, temperatures and air speeds. Over the years a great deal of such 
information has built up and has been proved in practice, so that there are 
accepted values of h for many purposes. An idea of the order of magnitude of h, 
and its variability with different materials, is given in Table 3. 


Table 3 
Process Type of boundary wim °C 
free solid—gas 0.5—1 000 
convection solid—liquid 100 —5000 
forced solid—gas 10 —1000 


convection solid—liquid 100 —10000 


Figure 2 


Aé is to be treated as a single 
symbol. It does not mean 
A x 0. 
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We can now go back to Figure 2 and write down an expression for the heat 
transfer right across from the warm air inside the room to the colder air outside. 
In the steady state, the rate of heat transfer will be the same at all points so that, if 
A is the effective area and b is the thickness of the wall, we can write 


q K 
A =h,,(6,, — 9,) = pala — 65) = hy (82 — Pout) 


where h,,, and h,,, are the convective heat transfer coefficients at the inside and 
outside surface respectively. 


Example 1 
By using the three heat transfer equations for Figure 2, show that 
ee eae 
= A(O.. —0,,.)4<—+-—+—> . 
q CA own tn 2 K + om 


Solution 
From the three equations, we have: 


q 

qb 
ae 
q 

02 — Gon = Fp 


Adding these equations, we get 


ee ee 
G.. ae a A ‘ie 56 K : —} 
Hence 
eee 2 
g= A(6.,, eed a) . K zs <j 


The expression which you were asked to derive in Example 1 gives the heat 
transfer rate in terms of the overall temperature drop. The quantity 


ke 
hin K Agut 


is usually referred to as the U value for the wall. The U value is the overall heat 
transfer coefficient for the wall and the surface effects combined. So the heat 
transfer rate can be written 


q ssa A U(6,, Be Cou) (2) 
where 
11 b Se ee, 
at 45 35) 
{~ = K - | 


In the technical literature, a = = +- = + _ is usually referred to as the total 


in out 


thermal resistance of the wall. 


Exercise 1 


A brick wall of a house, 0.24m thick, has a thermal conductivity of 0.7 Wm~!°C™~?. The air 
temperature on one side of the wall is 21°C and on the other side the temperature is 0° C. 
The convective heat transfer coefficients are 8 Wm~?°C~! for the high temperature side 
and 19 Wm~?°C~! for the low temperature side. Calculate the U value for the wall. 


Exercise 2 


For the wall specified in Exercise 1, estimate the heat transfer rate per square metre of wall. 
What is the thermal resistance of the wall? 


You met U values in the 
television programme. 
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Exercise 3 
If the resistance of the air film next to each side of the wall were to be neglected in 


1 1 
Exercises 1 and 2 (i.e. if it were to be assumed that "SE ok Cas 0), calculate the resulting 


out 


percentage error in the heat transfer rate. 


[Solutions to Exercises 1-3 on p. 28] 


Because of the varying area across which heat transfer takes place in pipes, U- 
values in the strict sense are not conveniently calculated, but an equivalent 
procedure is the subject of the following example. 


Example 2 
Derive an expression for the heat transfer rate in terms of the overall temperature 


difference for a pipe of circular cross-section carrying hot water and surrounded by 
cooler air. 


Sketch a graph showing the temperature variation between the hot water and the 
air. 


Solution 
I’ shall use the notation of Subsection 2.4 (ie. inner radius r,, outer radius r,, 
corresponding temperatures 6, and 0,). In addition I shall call the temperature of 
the hot water 0,,, that of the cool air 0,,,, and the convective heat transfer 
coefficients at the inside and outside surfaces h,, and h,,,, respectively. 
Then the equations for heat transfer by convection at the inner and outer surfaces 
give 

[= 2nr jth, Fin om 0,) igs 2arlh,,.(8, si Gia) 
while the equation for heat transfer by conduction through the pipe wall is given 
by Equation (5) of Subsection 2.4 as 


_ 2nkl(8, — 62) 


log, (r,/r,) 
Hence 
eee. ee 
ae 2nr, Ih,’ 
ee. Se 
ee! 2anr lho 
q log, (r2/r,) 
Adding, 


| 1 
ee + —log, (r9/r,) + ; 


ee r hin ae 


Thus, 


1 1 1 = 
q = 21(9i, — Pout) ic + “108. (r,/r,) + =| 


in 


The required sketch is Figure 3. 


4.2 Design for cooling 


Throughout most of this unit so far I have used examples in which heat transfer pie 


has, by implication, been presented as an undesirable phenomenon. Indeed the title 
of this section includes the word ‘insulation’, which simply means the cutting down 
of heat transfer rates. This is a very topical subject at a time when economy in the 
use of energy is, or should be, the order of the day, and I shall attempt to give it 
its due in the next subsection. In the meantime, however, it seems appropriate to 
remind ourselves that there are many occasions when heat transfer is desirable and 
ingenuity is properly exercised in fostering it. 
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One such example is an air-cooled internal combustion engine, such as most S 
motor-cycle engines. You must have noticed what the external cylinder wall of 
such an engine looks like—it has deep corrugations to increase the area of contact 
with the atmospheric air which acts as coolant. The same idea of an ‘extended 
surface’ is used in refrigerators and in industrial heat exchangers (which are 
heating or cooling devices) where, for example, pipes are fitted with fins as in 
Figure 4. 


I want to take a look at such a finned pipe, in order to model the effect 

of the fins. I shall assume that the pipe is carrying a hot liquid which is to be 
cooled by the air surrounding the tube. I propose to model a single fin on the 
assumption that its action is not interfered with by the other fins: in other words 
that each fin transfers energy to air at the ambient atmospheric temperature (not 
heated by the other fins). 


Figure 4 


Figure 5 


(ii) 

A single such fin is shown in Figure 5. Since its width is considerably greater than 
its thickness, I shall, for simplicity, neglect the heat transfer from its edges and 
model only the heat transfer from its top and bottom faces. This is a pessimistic 
assumption as far as the performance of the fin is concerned. We shall look at a 
small but finite length / of pipe and fin. Also, as before, we shall assume a steady 
state in which the temperature at any point does not change with time, and also 
that the temperature at any distance x from the pipe is uniform over the whole 
cross-section of the fin. 


We begin by considering the energy input to and output from a small element of 
the fin having width 6x at distance x from the pipe wall; that is to say the energy 
transferred to and from that element, which is shown in Figure S(i). Figure 5(i1) is 
a picture of the element in isolation. q, is the rate of energy transfer by conduction 
to the element from the part of the fin on the left of the element. q, is the rate of 
energy transfer by conduction from the element to the rest of the fin on the right 
of the element. q, is the rate of energy transfer (both sides) from the element to the 
atmosphere. Now if the temperature of the element is to remain constant, then 
during any interval of time, the amount of energy given to the element must equal 
the amount of energy removed from it. I mentioned this when I discussed the 
steady state in Subsection 2.2. Any net gain or loss of energy would result in a 
change of internal energy and hence in a temperature change. Thus we can write 


di: = 42 T Q3- (3) 


The next step is to try to express all the quantities in Equation (3) in terms of the 
dimensions and physical properties of the fin. 


For example, Fourier’s law tells us that 


where A = lb, 
k = thermal conductivity of fin material, 
0(x) = temperature of fin at distance x from pipe wall. 


It follows that 


dé 
Se Seal 


and 


dé 
q> —Kib- + Ox). 
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As for q3, I shall assume that all the values of 0 are low enough to make radiation 
inconsiderable and I shall treat it as the convective heat transfer rate from the 
element to the atmosphere. The area from which this takes place is 2/5x (top and 
bottom faces with the edge neglected). So if the atmospheric temperature is 6, and 
the relevant convective heat transfer coefficient is h, then 


q, ~h x 2ldx x (0 — 6,). (4) 


This equation is approximate because we have neglected (among other things) the 
variation in 0(x) across the element. The relative error in this approximation 
reduces as 0x is made smaller. 


Putting these expressions for q,, q, and q; into the ‘input-output’ formula (3) 
(rearranged into the form q, = q, — q,) gives 


dé dé 
2h(6 — 0,)6x = cb) (x + 0x) as (op (5) 


Since 6x is small, we can approximate the expression in square brackets by using 
dé 
the Taylor (tangent) approximation about x for the function = This 
approximation is 
dé dé d70 
ae + 0x) = a + ox a(x) ce 


and Equation (5) therefore simplifies to 
2 


d“@ 
21h(6 — 0,)dx =~ Klbox 5 See 


Dividing by 16x gives 
d?0 
AG = 0) ee (6) 
dx 


The sign ‘~’ and the ‘+...’ refer to approximations which become better and 
better the smaller 6x becomes. We can therefore make Equation (6) as accurate as 
we please by making 6x small enough. Taking the limit of vanishingly small 6x, we 
may therefore dispense with the approximation signs and write 


d*6 
2h(0 =< 6.) = Kb (7) 


although of course even this equation is not an exact model of the fin problem 
because of other approximations, already discussed, which do not depend on the 
size of Ox. 


Equation (7) is a linear second-order differential equation and can be solved using 
methods given earlier in this course. For conciseness we first define 


2h 
A =< as 
so that Equation (7) becomes 
d70 F . 
ae d= —A’6.. (8) 
Question 


What is the general solution of Equation (8)? 


Answer 


The general solution of 
5 ‘ 
ee A°O = —A’*O, 


is (as we saw in Unit 6) the sum of two functions: 


M101. MS283. 
Preparatory booklet for 
TM 281 students. 


See Unit 6 
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(i) any one particular solution, and 


(ii) the complementary function, which is the general solution of 
d’0 
—, — 70 =0. 
dx? 
Hence the complementary function is 
Fe + Ge *, 


where F and G are constants. It is easy to see that a particular solution of (8) is 


G = @.. 
Hence the general solution of (8) is 

0 = Fe* + Ge + 6, (9) 
Question 


Briefly, how do we find F, G and 0,? 


Answer 


6, has already been defined as the atmospheric temperature, which can be found 
by measurement. To find F and G we need two boundary conditions, which when 
substituted in (9) will give two simultaneous equations to solve for F and G. 


One of the two boundary conditions we need arises directly from the way in which 
we derived Equation (3), from Figure S(ii). Clearly, when x = w, q, = 0 because 
there is no part of the fin beyond x = w which could receive energy by conduction. 
Hence when x = w, 


and since x and A are non-zero it follows that = = 0 when x = w. 

Here, then, is one of our boundary conditions. The most convenient way to obtain 
the other one is to specify the temperature of the outside of the pipe wall, where 
the fin is attached. This can be measured without much difficulty. We are dealing 
with the steady state and so this temperature, like the others, will be independent 
of time. We therefore state our second boundary condition as 


0= 0, when x = 0 (0. constant). 


Exercise 4 


By using the boundary conditions stated above, derive F and G in Equation (9) in terms of 
0, 2 and w, and show that 


0— 0 


rice Ax (2 Aw — Ax) 
= —.— je" +e 


where T= 0, — 0,. 


Exercise 5 

The following data apply to a finned copper tube: 
w = 50mm h=10 Wm™~*°C"! 
b=2mm 9, = 150°C 
can Wa *C* 6. = 20°C. 


(i) Estimate the temperature midway along the fin. 


(11) Treating the temperature midway along the fin as the average temperature over the 
whole length of the fin, estimate the heat transfer rate from the whole fin per metre of axial 
length. 


[Solutions to Exercises 4 and 5 on p. 28 | 
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4.3 Insulation 


By insulation we mean a layer (or several layers) of material which does not easily 
transfer energy and which is deliberately introduced in order to reduce the rate of 
heat transfer. Very often insulation is used to save fuel. Lagging of pipes is one 
example: pipes are lagged by wrapping them around with layers of a material of 
low thermal conductivity. Still air is a very good insulator; this is why ‘cavity 
walls’ are successful in house-building: they consist of two layers of bricks 
separated by an air space. Double glazing with its air space between two glass 
panels works on much the same principle. Perhaps the commonest and best 
known example of insulation is provided by the clothes we wear. Their main 
function, at least in the British climate, is to maintain a layer of still air between 
the skin and the atmosphere in order to keep the skin temperature at a 
comfortable level. In very hot countries white clothes may also act as reflectors of 
radiation. 


However, things are not quite as simple as that. It is true that still air is a poor 
conductor, but as we have seen, air when heated may not remain still and is then 
capable of heat transfer by convection. This is why the insulating properties of, for 
example, double glazing do not necessarily improve as the thickness of the air gap 
is increased, as would be the case if conduction were the only heat transfer 
mechanism. Instead the variation of heat transfer by conduction and convection 
combined is as shown in Figure 6 in terms of a combined heat transfer coefficient 
h,. This one coefficient takes into account both conduction and convection in the 


, 1 
air gap. The resistance of the air gap can therefore be expressed in the one term : 


instead of the three that would be necessary if conduction and convection had to 
be treated separately. Clearly, as the air space is made thicker there is an 
improvement up to a point, after which it actually gets worse because of the effect 
of convection. The optimum thickness depends on the temperature difference 0 12 24 36 48 60 72 84 
across the air—for a temperature drop of about 33°C, for example, the best Air space thickness/mm 
thickness for the air gap is about 15mm. (I am assuming that the double glazing is 
fitted in order to provide thermal insulation only—its use for sound insulation is 
beyond the scope of this unit.) This effect, of convection preventing us from getting 
the full benefit of reduced conductivity, is found also in cavity walls. This is the 
reason why nowadays the cavity is often filled with plastic foam or granules. In the 
foam, and between the granules, there are many small air spaces, separated from 
one another so that there is no convection from one to the other. 


Figure 6 


The lagging of a pipe can also have its snags. One effect of winding on the 
insulating material is to increase the surface area in contact with the surroundings, 
and this tends to increase the heat transfer rate due to convection, although of 
course the surface temperature is reduced by the insulation. This point is taken up 
again in Problem 3 of Section 5. 


The following calculations are intended to illustrate some of the points I have 
made in this subsection, particularly the use of the combined 
convection—conduction coefficient. 


Example x. 


Assuming that.heating a house costs 16p per therm (1 therm = 1.055 x 10° J), 
work out the cost-of the energy which passes through 1 m? of window in one hour 
when the inside temperature is 18°C and the outside temperature 0° C: 


(i) for a single glass pane, 
(11) for double glazing. 
Assume the following data. 


Measurement data: glass panes are 6mm thick, with a 12mm air space for double 
glazing. 
Heat transfer data: « =1 for glass, units Wm~!°C"! 

h = 10 for the inside face, 


: —-20f7°:-1 
20 for the outside ae were 
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Solution 


(i) For a single glass pane (Figure 7), the energy transmitted per second per m7 
is given by three equations: 


q = 10(18 — 0,), 1e. 18 — 0, = = (inside convective heat transfer); 
i. 


0, -—@ 
q=1 | : 00 2, i.e. 0, — 8, = 0.006q (heat conduction through glass); 


q = 20(0, — 0), ie. 8, = = (outside convective heat transfer). 


Figure 7 
Adding these three equations: 


1 1 
18=q re + 0.006 + rH = 0.1564, 


10 20 
so that 


oe 
a 6156 


Thus in one hour, the energy transmitted through 1 m? of window is 


x 3600J = 415385 J. 
0.156 


10-8 
a ——— = (),0630 
Cost = 16 x 415385 x 1055 pence 
(ii) For double glazing (Figure 8), the energy transmitted per second per m7 is 
given by five equations: 


q = 10(18 — @,), ie. 18 — 0, = = (inside convective heat transfer); 


Figure 8 


gd, —8@ 
g=4 | ; 00 2, ie. 8, — 6, = 0.006q (heat conduction through inner 


pane); 


q =h,(0, — 3), ie. 0, — 03 = = (heat transfer through air space); 


0, a 0, * e 
q=1 0006 |}? ee 6, — 0, = 0.006q (heat conduction through outer 


pane); 
q = 20(0, — 0), ie. 0, = S (outside convective heat transfer): 


Adding: 


10 h 20 


c 


1 i 1 
18 = a[,5 + 0006 + + 0006 + 5) 


1 1 2 
q = 18 [7 + 0.006 + 7 +006 + 5] ; 


c 


In order to evaluate this we need to find the value of h, from Figure 6. This, in 
turn, requires a knowledge of the temperature drop in the airspace, i.e. 0, — 03. 
We can get an estimate of this by using a process of successive approximations. 


As a first approximation assume that the temperature difference across the cavity 
is (say) 6° C. Then, from Figure 6, for a 12mm airspace h, ~ 1.6. Substituting this 
into the expression for q we get 


q = 18(0.1 + 0.006 + 0.625 + 0.006 + 0.05)~? 


18 
= = 2G. 
0.787 : 
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Thus from our equation for 0, — 6, we have 


7 $6 


Looking again at Figure 6, the value of h, for a temperature drop of 14.3°C is 
about 1.7. This gives 


q = 18(0.1 + 0.006 + 0.588 + 0.006 + 0.05)~1 


18 
= —— = 24, 
0.750 he 
and 
24 
6, — 0, = — = 14.1. 
ee 


So, within the accuracy attainable with Figure 6, the answer is that there is a heat 
transfer rate of about 24 W through each square metre of window, and a 
temperature drop of 14°C across the air gap. (We are only looking for a 
reasonable estimate.) 


Thus in one hour, energy transmitted through 1 m? of glass ~ 24 x 3600J. 

Cost ~ 16 x 24 x 3600 ue 

ost ~ 16 x 24 x x 
1.055 

So double glazing has reduced the cost by a factor of 44 or so. Notice again that 
the calculation of the heat transfer rate across the cavity took into account both 
conduction and convection, by way of the combined convection—conduction 
coefficient obtained from Figure 6. The same applies to the exercise below, for 
which the combined convection—conduction coefficient is part of the data. 


= 0.013 pence. 


Exercise 6 


A brick wall of a house consists of two layers of brick separated by an air space. Work out 
the U value for this cavity wall and compare your answer with that to Exercise 1 on p. 18. 
Assume the following data. 


Measurement data: each layer of brick is 0.12 m thick. 

Heat transfer data: 
kK = (0.7 for brick, units Wm~/°C"! 
h= 8 for the inside face, 


= 19 for the outside s| units Wm~?°C™? 
h, = 1.6 for the cavity. 


[Solution on p. 28 | 


Summary of Section 4 


1. Convective heat transfer in a fluid can be modelled using the following 
assumptions: 


(i) the temperature in the bulk fluid is uniform, say 6,; 


(ii) if the fluid meets a plane solid surface, with area A and temperature 0,, then 
the rate of heat transfer from surface to fluid is hA(0, — 0,), where h is the 
convective heat transfer coefficient. 


2. The equation for heat transfer through a wall or window with area A is 
q = AU(6,, — 9.) 


where 6, and 6,,, are the bulk air temperatures inside and outside the building 
and U is a constant known as the U value. U values can be calculated by adding 
up the formulae for the differences between adjacent temperatures as in the 
solution to Example 1. 


3. Convective heat transfer from a surface such as a pipe wall can be increased 
by means of fins. The heat transfer in the fin can be modelled by assuming that 
the temperature in the fin depends only on its distance x from the pipe wall, and 
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then considering convection from each element of the fin and conduction from one 
element to another. This leads to the differential equation 

d70 2h 

= = (0 — 0) 

ax 0 


where h is the convective heat transfer coefficient, x is the thermal conductivity of 
the material of the fin, b is its thickness and 0, is the temperature of the air. 


4. The heat transfer properties of a layer of fluid (such as the air space in double 
glazing) are complicated, and best modelled by a combined heat transfer coefficient 
h,, which depends both on the thickness of the layer and on the temperature 
difference across it (see Figure 6). 


5 End of unit problems 


Read all the following five problems and attempt as many as you have time for 
without first looking at the solutions. Then compare your attempts with the 
solutions. Then read again the problems you have not attempted and read 
carefully through their solutions. In this way you may learn something from all the 
problems. (Solutions on pages 29-30.) 


Problem 1 
2kg of water at 80°C are added to Skg of water at 20°C and allowed to mix thoroughly. 


Neglecting any energy exchanged with the container or with the surroundings, estimate the 
final temperature of the water. 


Problem 2 


A storage vessel may be modelled as a hollow sphere of uniform material. The internal and 
external radii are r, and r, respectively. The temperatures of the inner and outer spherical 
surfaces are 0, and 6, respectively and the thermal conductivity of the material of the vessel 
IS K. 


(i) Show that the steady-state heat transfer rate through the wall of the vessel is given by 


‘eS ee 

ee 
(ii) Show that if the storage vessel is made of steel, and 0, = 20°C, 6, = 15°C and 
r, = 0.1m, then without insulation the heat transfer rate through the wall of the vessel 
cannot be kept below about 340 W however large the outside radius r, is made. 


(Hint: the surface area of a sphere of radius r is 4zr’.) 


Problem 3 


In the solution to Example 2 of Subsection 4.1 we found that the rate of heat transfer from 
the hot water inside an unlagged pipe to the cooler air outside is 


1 1 ee Be 
= 2n1(0,, — 95n)4—— + —1 ———> . 
q Tt ( in a + + K Og, (r,/r;) = =| 
Derive the corresponding expression for the heat transfer rate q,,, from a pipe fitted with a 
layer of lagging which has a convective heat transfer coefficient h,,, and coefficient of 
thermal conductivity x,, the new outside radius being r,. Hence derive the condition which 
the value of h,,, must satisfy if the new heat transfer rate is to be less than the original one. 


Problem 4 


This problem concerns the content of Subsection 4.2. By integrating the right-hand side of 
Equation (4) on p. 21 and using the result of Exercise 4 on p. 22, estimate the total heat 
transfer rate per metre of axial length from the fin specified in Exercise 5 to the atmosphere. 
Compare your result with the solution to Exercise 5(ii) and comment on the difference. 


Problem 5 


The average rate at which energy is received from the sun at the surface of the earth is 
about 1347 Wm”. Estimate the surface temperature of the sun, assuming (i) that the 
constant p in the radiation equation is proportional to the surface area A of the radiating 


Energy 


af = ; ; a a > radiated from 
body, (ii) that for the sun uy = A x 56.7 x 107? (the units being W K ~*), (iii) that the ratio this part of 
of the radius of the earth’s orbit to the sun’s radius is 216, and (iv) that the earth’s radius the surface 


may be neglected. 


(Hint: the beams of radiation from the sun’s surface spread out as indicated in Figure 1.) Figure 1 
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Appendix 1: Solutions to the exercises 


Solutions to the exercises in Section 1 


1. From Table 1, c = 4200, and from Equation (1), 
E, — E, = mc(6, — 6,) 
= 2 x 4200 x (80 — 20) 
= 504000 
= 5.04 x 10°, 
i.e. there is an increase in internal energy of 5.04 x 10°J. 


2. The energy output of the element is 2000 Js~!. Thus the 
time to produce 5.04 x 10°J will be 


5.04 x 10° 
2 x 10° 

= 252 seconds 

= 4 minutes 12 seconds. 


seconds 


3. Some of the energy is used to heat the kettle, and some 
is used to heat the air surrounding the kettle. The rate at 
which energy is used to heat the water is therefore less than 
was assumed in Exercise 2, and so the time obtained is an 
underestimate. 


—--—————————————————————— 


Solutions to the exercises in Section 2 


1. From Table 2, x = 0.7. Substituting this and the values 
A = 1, b = 0.24 into Equation (2), we obtain the rate of heat 
transfer through 1 m? of wall: 


18 — 10 
2654451 
ae «| 0.24 | 


== 23,3, 
Thus the rate of heat transfer through 1 m? of wall is 23.3 W. 


—8 
2. The temperature gradient is a — 33.3, the units 


being °Cm~*. Let the temperature at 100mm from the inner 
face be °C. Then 


6 — 18 
eS ae 
0.1 ie 
9 =18 — (01 x 33.3) 


~ 14.7, 
and the required temperature is 14.7° C. 


3. From Table 2, x = 1.0 for glass. Substituting this and the 
values for |, 6,, 6,, r; and r, given in the exercise into 
Equation (5), we obtain the rate of heat transfer through one 
metre length of pipe wall: 

_ 2x x 1.0 x 1 x (90 — 50) 

= log, (40/25) 

= 535. 


Thus the rate of heat transfer through 1 m length of pipe wall 
is 535 W. 


4. Let the temperature be 0°C at radius x = 0.032. Then by 
Equation (4), 
a3 


Be ee ane 
SS CE aes ne aa 


and we can find C because we know that 0 = 90 at 
x =r, = 0.025,. 


Thus, 
535 
= ———______] i 
= 2x x 1.0 x 1 eo © 
= 314+, 
so C= —224, giving 
0 = 69.1, 


ie. the temperature in the pipe wall at a radius of 32mm is 
69.1°C. 


An alternative method would be to observe that Equation 
(5) applies to that part of the pipe wall lying between 

r, = 25mm and x = 32mm, so that 

2x x 1.0 x 1 x (90 — 8) 


aes log, (32/25) 


een er SSS SSS se 


Solutions to the exercises in Section 3 


1. In this case, the length of the bar (0.2m) is the quantity 
corresponding to the thickness b in Equation (2) on p. 11 of 
Subsection 2.3. The cross-section area in m7? is 


7 
A =— x 0.012? 
4 


= 413i = 10-*. 
The rate of heat flow in watts is g = 1.12. 
Since 6, — 0, = 10, Equation (2) gives 
Wes a 
giving 
kK = 198, 


i.e. the thermal conductivity of the material is 198 Wm /°C~'. 


2 @ -~250°C = 073 — 250K = 23K. 
Gi) 3°C = (073 + 3)K = 276K. 
3. For radiation, q = 64, = u(0 + 273)*. 


Assuming that the filament dissipates its energy by radiation 
only, 


40 = n(0 + 273)*. 
If 8 = 2000, then 

40 40 
~ (2273)* 2.67 x 1083 
the appropriate units being W K ~¢. 


u = 1.5 x 10712, 
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Solutions to the exercises in Section 4 


7 b aa 
I U=4—+-+— 
K 


hin hous 

‘ ‘5 ee at 
Be 07: 19 

= 1.92, 


the units being the same as those of h, namely Wm~?°C™?. 


2. For the rate of heat transfer through 1 m? of wall we put 
A = 1 in Equation (2), and obtain 


q=1 x 1.92 x (21 — 0) = 40.3. 
Thus the rate of heat transfer through 1 m? of wall is 40.3 W. 


oe 
The thermal resistance is — = = 0.521, the units here 


1 
See BP 
being m?°C W~?. 


i 1 ; 
3. Putting "es Ses 0, we obtain 


mm out 
2 = 
Ag? 
= 2.92. 


The percentage error im the heat transfer rate will be the 
same as the percentage error in U, namely 


2.92 — 1.92 


Top X 100 = 52.1% 


4. By Equation (9) on p. 22, 
0 — 6, = Fe* + Ge*. 
At x = 0, 6 = @, (one of the boundary conditions), 


SO 

0, -6,=F+G=T. (1) 
Also, 

“(0 — 0,)= . = AFe — 1Ge~*. 
Atx=w, = = 0 (this is the other boundary condition). 
Hence, 

AFe*” — AGe"*” = 0. (2) 
From (2), 

fg sie 
Le. 


Fo™ = G. (3) 
Substituting this expression for G into (1), |: 
Fi +e™=T. 


Thus, F = 


and using (3) again, 
Te?2*” 
= 4+ e24w 


Hence, 0 — 6, = ie? pee) 


ie 


as required. 
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5. 


[2h 
ge ee 
(i) Kb 
= 240 ae 
Al 380-x 0002 °° 


We also have T = 150 — 20 = 130, w = 0.05, and midway 
along the fin x has the value 0.025. Hence, 
130 
g. = ———_____ x 
* 1+ exp(2 x 5.13 x 0.05) 

x {exp (5.13 x 0:025) + exp (2 x 5.13 x 0.05 — 5.13 x 0.025)} 
ae 
~ 1+ 1.670 
== 327. 


Since 0, = 20, we conclude that the temperature half-way 
along the fin is 147°C. 


(11) Assuming 147°C to be the mean temperature of the 
fin, the heat transfer rate is given by 


gq = 2whh(Onean — 94) 
(where the factor 2 occurs because the fin has two sides). To 


find: the heat transfer per metre of axial length we put / = 1. 
Thus: 


gq =2 x 0.05 x 1 x 10 x (147 — 20) 
= 127, 


so that the heat transfer rate from the fin is 127 W per metre 
length. 


a 


{1.137 + 1.469} 


6. A section of the cavity wall is shown below. 


Inside 21 °C Outside 0°C 


pss] 


Using the notation of this figure, the energy transfer per 
second per m7? is given by five equations: 


for the inside layer of brick, 


g= 8(21 _ 0.) i.e. 21 — 6. — = 
| See 0.12¢ 
— 0.7 baba 2 i. 6. — = ——: 
q | FE , ie. 6, — 0, 07” 
for the cavity, 
q 
= + Q — Q sees = ee 
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for the outside layer of brick, 


0, 5 : 0.12q - 
q ao 0.7 | 0.12 ; L€. 0, Sa a... = 0.7 9 
q 
q= 19/0 out oe 0), 1.¢. is 19° 
Adding 


21 = Mien ce oe 
Sle Oa “46° nr er 


We oo ee 
= 0.873. 
That is, the U value is 0.873 Wm~?°C7?. 


-1 
The pe f 5 O12. 1 OL 1 
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In Exercise 1, for the same total thickness of brick but 
without the cavity we have a U value of 1.92 Wm~?°C™/, 
so the cavity has reduced the heat transfer rate to less than 
half of the former value. Filling the cavity with plastic foam 
can reduce the U value to about 0.5 Wm~?°C™?. 


Appendix 2: Solutions to the problems 


1. Heat transfer will cease when all the water is at the same 
temperature, say 6°C. 


Then the change in the internal energy of the warmer water 
is {2 x 4200 x (0 — 80)}J, and the change in the internal 
energy of the cooler water is {5 x 4200 x (@ — 20)}J. Since 
no energy is given to or received from the surroundings, the 
total internal energy of the water must remain constant, and 
hence the total change in internal energy must be zero. Thus, 


{2 x 4200 x (0 — 80)} + {5 x 4200 x (@ — 20)} = 0. 


Cancelling out the factor of 4200 and collecting terms, we get 


70 = 260, 
6 = 489 — 37.1. 
Thus the final temperature will be 37.1°C. 


2. 


(i) The temperature at any point in the wall is determined 
solely by its distance from the centre of the sphere. Hence the 


_ ad 
temperature gradient in the steady state is gs where @ is the 
x 


temperature at a point whose distance from the centre of the 
sphere is x. The rate of heat transfer through the entire 
spherical area at distance x from the centre is given by 


gq = —K Xx 4nx?” x —. 


D | 


Figure 1 


Hence 
a0. 
dx  4nxx?’ 


In the steady state, g is a constant, so this equation can be 
solved by direct integration to give the general solution 


—q 
= d 
‘ | Po os 


ae 
4nKx 


where C is a constant of integration. 


Thus since 9 = 6, when x = r, and 6 = 0, when x = r,, we 
have 
se 


Yo MERE, 


=e 
ArKr, 


2 


and subtracting, 


so that 
1 i ae 
q = 4nk (8, — 0) _ | 
i pee 


(ii) As r, increases, the heat transfer rate derived above 
decreases, to a limiting value of 


1 eh | 
fel 
For the given figures (and taking x for steel = 54 from Table 


20n p. 10), this gives a heat transfer rate of 339 W. Thus the 
heat transfer rate cannot be kept below about 340 W. 


3. A section of the lagged pipe is sketched below. 


Figure 2 Bout 


Lagging 


We have four equations for heat transfer from a length / of 
this pipe: 


: ia 
Gisg = 207th, (6,,-—9,), ie. 6,,—0, = ro The 
2nkl(6, — 9) : Aiag 108, (T/T) 
= —_____— e. 6, —0, = ss 
iag log.(r2/r1) = : 2nkl 
eS 2k (0, — 93) — Nag log, (73/r2) | 
diag = log.(r3/r>) ee 2nk,l 
: ia 
Niag = 2nrslh,,.(83 = 03%, 1.€. 6; — a —— te 
Adding, 
5° oe 1 “ log, (r,/T;) x log, (r3/r2) a 1 
a K Ky Flin.) 
and so 


diag = 271(6., — Pour) X 


r,h., K Ky r3h 


| 1 _ ei) OR 1 . 
lag 


Thus, to meet the condition q,,, < q we must have 


1 4 Rat) OR) 1 


PM K Ky r3Mae 
1 es log, (r,/r;) re 1 
ry hin K rhea 
that is, 
l 1 1 
ont) VA 
Ky r 3 Mag rohgut 
t 1 eee). 
r3Mag rohgut K 1 


1} : ee Cala)} 
3 ohour Ky 


This means that if h,,, is too high the heat transfer rate of 
the lagged pipe can, at least in principle, be greater than that 
of the unlagged pipe! 
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4. By Equation (4) on p. 21, the rate of heat transfer from 
an element of the fin of width 6x is given by 


h x 21l6x x (6 — 6,). 
Thus the total heat transfer rate for the whole fin is 


| 2Ih(6 — 6,) dx 
0 
“ 2h | (0 — 6.) dx 
0 


ae 2 
= 2Ih | ——- {e* + 4") dx (by Exercise 4) 
o ite 


2IhT bs ’ ge I 


eS eee 
Cs ae 
~ qa em te 
_ 2hT (e*” — 1) 
A(e?4 + 1) 


0 


isis e*”) ins (1 es "pices | 


For the figures in Exercise 5 this gives, for ] = 1: 
2 x 10 x 130 x (1.670 —1) | 127 
5.130 x (1.670+1) 


giving a heat transfer rate of 127 W per metre length. This 
agrees (to three figures) with the result of Exercise 5. 


5. Below is a sketch (not to scale) of a beam of energy 
radiated from the sun to the earth. 


Figure 3 


If the beam has a cross-sectional area of 1m? at the earth’s 
surface, it will have a cross-sectional area of (34<)? m? at the 
surface of the sun. Thus the rate of radiation from (54)? m? 
of the sun’s surface is 1347 W. Hence, 


1347 = p x (6 + 273) 
where pt = (s4<)? x 56.7 x 10-9. This gives 
1347 x (216)* x 10° 
5.67 
= 11083886 x 10°. 


(0+ 273y = 


Thus, 
6 + 273 = 57.7 x 10* = 5770, 


giving a surface temperature of 5497°C (or approximately 
5500° C). This is, of course, only a rough estimate. 
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